Specific heat at tlie transition in a superconductor with fluctuating magnetic moments 
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In the heavy-fermion materials CeCoIns and UBeia, the superconducting order parameter is 
, coupled to flucutating magnetization of the uncompensated part of the localized /-moments. We 

find that this coupling decreases the superconducting transition temperature and increases the 
jump of the specific-heat coefficient, which indicates entropy transfer from the magnetic to the 
, superconducting degree of freedom at the transition temperature. Below the transition, we find 

D . that the magnetic fluctuations are suppressed. We discuss the relation of our results to experiments 

pLn ' on CeCoIns under pressure. 
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In this paper, we examine the effects of fluctuating paramagnetic moments on the specific-heat jump at the transition 
^ AC. We are motivated to examine this question due to recent experiments in CeCoIns, a heavy-fermion supercon- 
I , ductor with an unusually large ratio of AC/jTc « 5 at ambient pressure^, whereas the BCS value for singlet s-wave 
^ ' superconductor is 1.43. Anisotropy of the gap further decreases this value; for instance, for a d-wave or an anisotropic 
^ p-wave superconductor, it is suppressed by factors 2/3 or 5/6, respectively. We also recall that UBeia^, one of the 
first heavy Fermion superconductors to be discovered, has AC/jTc « 2. In UBcia, the normal-state specific- heat 
coefficient 7 above Tc increases as temperature is lowered indicating that the asymptototic Fermi-liquid state is not 
reached down to Tc- In CeCoIns, the integrated measured entropy 
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C/TdT (1) 



is about 3 times larger than '^{Tc)T^ assuring that were it not for the superconducting instability, a further enhance- 
ment of the normal state specific heat at lower temperatures would occur. Indeed, in a magnetic field just above 
K>" ' Hc2 = 5T, the specific-heat coefficient has a sharp upturn around IK as T decreases below Tc and continues to rise 
i down to the lowest measured temperature 0.2K. The entropy is quite adequately balanced between this measurement 
^2 ■ and the zero-field entropy of the superconducting state demonstrating the consistency of the experimental measure- 
' ments. We may associate the high temperature constant 7 state to an unstable Fermi-liquid state characterized by 
• a Fermi temperature Tph ~ 35 — 50K (see (I1U|I below), and the low temperature constant 7 state to an asymptot- 
fT^ ' ically stable Fermi-liquid state characterized by a temperature Tpi- The measurement of 7 above 5T suggests that 
. Tpi ~ IK. Thus, Tph ^ Tc- We see that Tc is bigger than Tpi by a factor 2 or 3. To simplify the analysis, we shall 
assume Tc ^ Tpi. We shall see that the results still give a good agreement with the experimental data. 

The two energy scales can originate in heavy fermions from a variety of causes. For example, multistage Kondo 
effects can arise in magnetic impurities when crystal field splittings are taken into account^; the values of the crystal- 
field splittings have been measured in various experiments^'^. In fact, a low Kondo temperature scale of about 1.5K 
has been deduced^ from the crystal-field scheme and the high Kondo temperature of 35K. A suppressed putative 
Q ] antiferromagnetic transition below Tc could also provide a lower energy scale. Such a transition would of course also 
Q . be suppressed in a field as large as 5 Tesla required to suppress superconductivity and the j{T) would then show a rise 
■ as observed below about 1 K (followed by a decline). In a recent detailed analysis®, the specific- heat and magnetic- 
. ^ , susceptibility data of CeCoIns have been decomposed into a contribution from a heavy Fermi liquid characterized by 
' Tph and a contribution from local impurities with Kondo temperature Tpi . The microscopic origin of the scale Tpi 
JH \. is not important in our model and the ensuing calculations. 
. . . . We shall study the renormalizations of the specific heat (and magnetic fluctuations) using the following Ginzburg- 
Landau free energy function. 

^[^,M(r)] = I d'^rl^^^^l^f +^-\^\'+ "^^ + V ^(r) + ||^|^M^(r) 
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^^V r + 7lV^l1 + y^HT + T,f) + vmUi. (2) 



In (|21), i/" is the superconductivity order parameter, M(r) is the magnetization at r and the coupling rj between the 
superconducting order parameter tp and the magnetization M(r) is of the lowest allowed order and is independent of 
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momentum and frequency.The magnetic fluctuations are assumed to be local and classical: 
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(3) 



The assumptions underlying the free energy JSJ are the following: The fluctuations characteristic of the high tem- 
perature renormalized Fermi energy Tph are assumed to lead to pairing and to determine Tc and the bare parameters 
a and b in Eq. The fluctuations characterestic of the low temperature Fermi energy Tpi have a characterestic 

energy Tgf = 0{Tpi). As stated above, we further assume Tph ^ Tpi for CeCoIus. Two simplifications follow 

from this which are used in (jSJ: (i) The magnetic fluctuations characterised by Tgf can be treated classically near Tc, 
(ii) there is clear separation in the microscopic processes leading to T^, a, b on the one hand and those parametrized 
in Tj. Note that the term involving 77 correctly describes the Abrikosov-Gorko\ii effects of magnetic impurities on 
superconductivity to leading order. 

We can obtain the renormalized superconducting transition temperature, the jump of the specific-heat coefficient, 
and the susceptibility from the Helmholtz free energy. To calculate the Helmholtz free energy, we must integrate 
over all possible configurations of M and ^jj. We shall do this in a two-step process. First we integrate over M, 
because the integral is Gaussian, and can therefore be done exactly. That gives us a contribution .F/juctdV'l^) to the 
superconducting Ginzburg-Landau free energy. We then approximate the integral over tj) by the value of the integrand 
at its maximum. 
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is the free energy of the uncoupled magnetic fluctuations, and / is the lattice constant. This formula is correct in 
the classical regime of the magnetic fluctuations, Tpi << T that we are studying here. In the quantum regime, that 
is, when the temperature is lower than the energy of the magnetic mode at momentum kc, we have to cut off the 
momentum integration at the value corresponding to T. 

The Ginzburg-Landau energy density for the superconducting order parameter after including the effect of the 
fluctuating magnetization is 
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Here we have combined 77 and a into one parameter with the dimension of temperature 
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which characterizes the strength of the magnetic fluctuations. From the functional l(21), we see that 
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is the jump of the specific heat coefficient at Tc of the bare superconductor per unit cell, so it has dimension of inverse 
temperature. Its value T-y is of the order Tph- From the experimental values of 7(T > Tc) ~ 0.2 — 0.3J/molK^ using 
R ~ 8J/molK, we obtain 



~ Tph = 35- 50K. 



(10) 
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FIG. 1: (a) Renormalized transition temperature and (b) renormalized jump of the specific-heat coefficient as a function of 

Once T-y is more than an order of magnitude bigger than the transition temperature, the measurable quantities are not 
very sensitive to its precise value. Behavior of F^c close to zero defines the renormalized Ginsburg-Landau parameters 
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The new transition temperature T* satisfies the equation 



6n = t: -t, = -t, 



and the new specific heat jump is given by 
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We note that the coupling to the magnetic fluctuations decreases the transition temperature and increases the 
specific-heat jump. We would like to compare our results to the experiments on CeCoIn5 under pressure'^, which show, 
that applying pressure increases and decreases the jump of the specific-heat coefficient. In the first approximation, 
we will assume that the bare superconducting parameters a, 6, and Tcj ^ well as the bare magnetic parameters a 
and Tsf are set at high energies, and thus do not change with pressure. The pressure is expected to change only the 
low-energy coupling 77 and hence the temperature T^m- In Fig^ we plot the dimensionless quantities T*/Tc and 
A7*/A7 as functions of the dimensionless variable T^m/Tc, to show that T* decreases and A7* increases as T^m 
increases. Hence, our interpretation of the experimental results is that the hydrostatic pressure suppresses T^m- In 
Fig^ we used the values Tc = 6K, T,/ = 1.5K (in agreement with the measurements of 7 in the magnetic fieldi 
discussed in the introduction), and T-y = 50K. Note that T^f < Tc, which justifies our assumptions that the fluctuating 
moments are classical. These values are reasonable on physical grounds, and give a reasonable agreement between 
the experimental data and the calculated results H12|l and (|13l) . Experimentally, as the applied hydrostatic pressure 
changes from to 1.48 GPa, T* and A7*/A7 change from 2.3 to 2.7K and from 3 down to 2, respectively. The 
experimental data correspond to T^m ~ T^. From the formulas (|12|l and H13|) . we find that as T^m/Tc changes from 
1.05 down to 0.85, T* and Aj* /Aj change from 2.2 to 2.7K and from 3.0 down to 2.1, respectively. The fit to the 
experimental data is, of course, not exact because at each value of the pressure, we are solving two equations 111(1 
and H12|l for one unknown T^a/ . For a detailed comparison with the experiment, we also need to assume a pressure 
dependence of another bare parameter, for example Tsf, but that pressure dependence comes out much weaker than 
the pressure dependence of T^m, in agreement with our considerations at the beginning of this paragraph. 

Finally, we calculate the change of the magnetic susceptibility relative to its bare value © due to the interaction 
77 to see the effect of superconductivity on magnetic fluctuations. The magnetic susceptibility is given by 
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we do both integrals by the two-step process described earher. Above T*, -0 = 0, so the susceptibihty stays equal to 
©. Below T*, the M integral gives 

y^"-l"'l'' " - 2(a(r^r.,) + ,,N.p) n( „,T..,,H.,M- j ■ <i^) 

The extra factor 1/ {a{T ~\-Tsf) +?7(fc)|V'|^) on the right-hand side will change the extremum in the numerator of p4|) 
by a number of the order of 1/V relative to the denominator, so we will neglect that difference. Then everything but 
this factor cancels between the numerator and the denominator, and we get 

'^^ = a(r + T,/) + r7|^(T)|2 ^^^^ 

with \il}{T)\'^ obtained by minimization of We see that the magnetic fluctuations are suppressed by coupling to 
superconductivity, but susceptibility of the localized moments remains finite down to zero temperature. Both effects 
are consistent with measurements of Knight shift^. 

In conclusion, we propose that the anomalously big specific-heat jump at the superconducting transition in UBeia 
and CeCoIns is caused by interaction between superconducting ordering and fluctuating unordered classical magnetic 
moments of localized f-electrons. We describe such a coupling between the two degrees of freedom by the simplest pos- 
sible Ginsburg-Landau functional. We find that the transition temperature decreases and the jump of the specific-heat 
coefficient increases as the strength of the paramagnetic fluctuations and/or of their coupling to the superconducting 
order parameter increase. Below the transition, magnetic susceptibility is suppressed but remains finite down to zero 
temperature. 
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